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We study a mean-field model for a clustering process that may be 
described informally as follows. At each step a random integer k is 
r~^ I chosen with probability pk, and the smallest cluster merges with k 

^^ ■ randomly chosen clusters. We prove that the model determines a con- 

\l I tinuous dynamical system on the space of probability measures sup- 

f">». ■ ported in (0, oo), and we establish necessary and sufficient conditions 

^D I for approach to self-similar form. We also characterize eternal solu- 

QP ■ tions for this model via a Levy-Khintchine formula. The analysis is 

based on an explicit solution formula discovered by Gallay and Mielke, 
extended using a careful choice of time scale. 
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1 Introduction 



1.1 A mean-field model for clustering 

The clustering processes we consider are motivated by a simplified model 
for domain wall motion in the one-dimensional Allen-Cahn equation dfU = 
dxxU + u — u^. The domain walls become points on the line, and the domains 
are the intervals separated by these points. The pattern coarsens by a simple 
rule: At each step, the smallest domain combines with its two neighbors 
to form a single domain, and this is repeated indefinitely. Computational 
simulations of this 'min-driven' domain coarsening process indicate that for 
a considerable variety of initial distributions, the domain size distribution 
approaches self-similar form [4, 12]. 

A mean-field model of this process was derived by Nagai and Kawasaki 
[12], and it turns out to be amenable to a rigorous analysis aimed at ex- 
plaining this behavior [4, 8]. We consider an infinite number of domains 
on the line, and study the statistics of domain sizes using a number den- 
sity function f{t,x). We assume that in any interval / of unit length, the 
expected number of domains with lengths in the range (x, x + dx) is given 
by f{t,x)dx. The expected value of the total number of domains in / is 
denoted N{t) = L f{t,x)dx. We assume that N{t) is finite, and denote 
the associated probability density by 

P.W = ff . (1.1) 

We let l{t) denote the size of the smallest domain at time t, so f{t,x) = 
for X < l{t). The expected number of coalescence events per unit time is 
then 

f{t,i)i. (1.2) 

The coalescence events affecting domains of size x in the time interval {t,t + 
dt) are (a) loss: consecutive domains of size x, I, y (or y, I, x) combine to form 
a domain of size x + l + y; (b) gain: domains of size y,l,x — y — I combine to 
form a domain of size x. Under the mean- field assumption that coalescing 
domains have sizes chosen randomly and independently from the current 
overall size distribution, these events have respective relative probability 
density 

Pt{x)pt{y), pt{y)pt{x), pt{y)pt{x-y-l). 
The rate equation for the evolution of / is obtained by summing over all 
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loss and gain terms: 



dtf{t,x)= (1.3) 

/ j'X — l POO \ 

fitJ)ii Pt{y)pt{x-y -l)dy -2pt{x) l pt{y)dy], x > I. 



\Ji Ji / 

Clustering phenomena are seen in fields as varied as population genetics 
and physical chemistry. Thus, while coarsening of intervals provides concrete 
motivation, the notion of a 'cluster' can have widely different interpretations 
in applications. The model (1.3) above is one of a family of what we call 
min-driven clustering models that can be analyzed together in one setting 
as in [8]. At each step a random integer k > 1 is chosen with probability 
Pk, and the smallest cluster merges with k randomly chosen clusters. The 
mean-field assumption is that all these random variables are independent. 
The only assumptions we impose on the probabilities are that 

oo oo 

Pk>0, ^pfe = 1, '^kpk <oo. (1.4) 

fc=l k=l 

The evolution of the number density under this process is described by the 
following rate equation for cluster size density, obtained by summing over 
all loss and gain events as in (1.3): 

oo 

dtf{t,x) = f{t,l)iY,Pk{pf{x-l)-kpt{x)), x>l{t). (1.5) 

k=l 

Here the notation pf' denotes /c-fold self-convolution. The case pi = 1 
corresponds to a "paste-all" model discussed by Derrida et al. [6]. Also see 
[9] for a model of social conflict with rather similar solution formulas. 

An important feature of the model (1.5) is its invariance under reparame- 
trization in time. If we change variables via t = T{t), f{x,t) = f{x,t), 
l{t) = l{t), then equation (2.5) retains its form since 

d'j{x,i) = fdtf, dii = fdti. 

A careful choice of the time scale is key to the analysis. In the first math- 
ematical study of (1.3) [4], the authors imposed the relation f{t,l)l = 1, 
meaning the number of coalescence events per unit time is constant [4]. In 
a more recent paper, Gallay and Mielke parametrized time by the minimum 
size, that is l{t) = t [8]. This leads to an elegant solution procedure that 
was used to prove some basic results on well-posedness and the approach to 
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self-similarity. Gallay and Mielke showed that (1.5) defines a strongly con- 
tinuous flow in L^, that (1.5) admits a one-parameter family of self-similar 
solutions and that suitable initial densities yield convergence to self-similar 
form. Precise comparisons between these results and ours are made later in 
this paper. 

In this article, we introduce yet another time scale. We parametrize time 
inversely to the total number of domains, so that 

t = N{t)-^. (1.6) 

We shall argue that this is a natural choice for a number of reasons. It 
retains the simplicity of the choice l{t) = t, and allows us to obtain: (a) 
existence and uniqueness for measure- valued solutions; (b) necessary and 
sufficient conditions for convergence to self-similar form; (c) a characteri- 
zation of eternal solutions for the dynamical system defined by (1.3). We 
comment on these in greater depth below. 

Let us first explain one simple motivation for (1.6). We show below that 
for (1.3), fit, l)i = -N/2. Thus, /(t, 1)1 = N'^ /2 when (1.6) holds, and (1.3) 
now takes the form 



dtf{t,x)= (1.7) 

' j'OO 

fit,y)fit,x-y -l)dy - f{t,x) f{t,y)dy, x>l. 



-1 fX — l 



2 

If we take l{t) =0 (as a model when the smallest domains have negligi- 
ble size, for example), this reduces to a basic solvable model of clustering: 
Smoluchowski's coagulation equation with constant kernel. In recent work, 
we provided a comprehensive analysis of dynamic scaling in this equation 
by exploiting an analogy with the classical limit theorems of probability 
theory [10, 11]. We use these insights to guide our study of (1.5). 

1.2 Measure-valued solutions 

Mean-field models of domain coarsening, such as the LSW model, Smolu- 
chowski's coagulation equation, and (1.5), correspond to physical processes 
where mass is transported from small to large scales. For several reasons, 
it is natural to consider measure-valued solutions, for which the size distri- 
bution need not have a continuous or integrable density. Such solutions are 
physically meaningful, as many clustering processes (e.g., polymerization) 
involve a discrete set of sizes based on an elementary unit. A formulation via 
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measures is also mathematically elegant, as it allows us to unify the treat- 
ment of discrete and continuous coagulation models, exploit simple criteria 
for compactness and continuity, and prove basic uniform estimates. 

To any solution of (1.5) we associate a probability measure Ft with dis- 
tribution function written 

^t(^) = ^y^ f{t,y)dy = J^ Pt{y)dy. (1.8) 

For any probability distribution F on [0, oo), we call / = inf{a;lF(x) > 0} 
the min of F. (This is short for "minimum size," regarding F as a probabil- 
ity distribution for size.) We will prove that the initial- value problem for an 
appropriate weak form of (1.5) is well-posed for initial probability measures 
FtQ with positive min. That is, (1.5) with (1.6) determines a continuous 
dynamical system on the space of probability measures with positive min, 
equipped with the weak topology. See Theorem 4.3 below. By comparison, 
Gallay and Mielke established that the initial-value for (1.5) defines a con- 
tinuous dynamical system on the space of probability densities in -L^(l, oo) 
equipped with the strong topology [8, Thm. 3.3]. The solutions we construct 
arise by a natural completion of these L^ dynamics. 

1.3 Dynamic scaling 

A common theme in recent studies of dynamic scaling in mean-field models 
of coarsening is that the approach to self-similarity is both degenerate and 
delicate. The problem is degenerate because there is a one-parameter family 
of self-similar solutions. For the model studied here with pj^ = Q for all k 
large enough, Gallay and Mielke found a family of self-similar solutions that 
may be rewritten in the time scale (1.6) in the form 

Ft{x) = F'^') {j^^ , l('\t)=t'/', ^G(0,1], t>0. (1.9) 

Here F^ ' is a probability distribution with density p^- ' supported on [1, oo). 
The density p^^' is known explicitly only through its Laplace transform. 
Only p'^) has finite mass (first moment); p^^^'{x) decays exponentially as 
X -^ oo. The distributions F^- ' for < 9 < 1 have heavy tails, with 
p(^)(x) ~ C0X^^^^^' as X ^ oo (see Theorem 5.1). 

The problem is delicate because the domains of attraction of the self- 
similar solutions are determined by the tails of the initial size distribution, in 
the precise manner explained below. (See [10, 13] for analogous results on the 
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LSW model of Ostwald ripening and Snioluchowski's coagulation equations 
with solvable kernels.) Gallay and Mielke showed that all densities with 
finite mass are attracted to the self-similar solution with 6 = 1. Moreover, 
for < ^ < 1 they showed that if the initial data pt^ is sufficiently close to 
p^^' in a suitable weighted norm then the rescaled probability density tpt{tx) 
approaches p^ ' with a rate of convergence determined by the weighted norm 
(see [8, Thms. 5.5, 5.7]). These results provide sufficient conditions for 
approach to self-similarity. 

Our aim is to establish conditions that are both necessary and suffi- 
cient to answer a more general question about arbitrary scaling limits. We 
characterize the set of all non-degenerate limits under a general rescaling of 
the form Ft{X{t)x) where A(t) is a measurable, positive function such that 
limi^oo A(t) = cx). A limit is non-degenerate if it is suitable data for the 
initial- value problem. That is, non-degenerate limits are probability distri- 
butions with a positive min. 

Theorem 1.1. Let Iq > 0, let Ft^^ be an arbitrary probability measure on 
(0, cxd) with positive min, and let Ft (t > t^) be the associated measure- valued 
solution of (1.5) (see Theorem 4-3)- 

(i) Suppose there is a measurable, positive function A(t) — > oo as t ^f oo 
and a probability measure F^ with positive min, such that 

lim Ft(A(t)x) = F,(x) (1.10) 

at all points of continuity of F^,. Then there exists 6 € (0,1] and 
a function L slowly varying at infinity such that the initial data Ff^ 
satisfies 

/ yFt„{dy) ^ X ~ L{x), as x — > oo. (l-H) 

Moreover, the min l{t) of Ft and the rescaling X{t) satisfy 

\{t)h^l{t)^t^^^L{t), as t^oo, (1.12) 

where l^ > is the min of F^ and L is slowly varying at infinity, 
related to L by (6.9). 

(a) Conversely, assume there exists 9 G (0, 1] and a function L slowly 
varying at infinity such that the initial data satisfies (1.11). Then l{t) 
satisfies (1.12), and 

lim Ft{l{t)x) = F'-^^x), X G (0, oo). (1.13) 

t— >oo 



Dynamics and self-similarity in min-driven clustering 7 

A positive function L is slowly varying at infinity if it is asymptotically 
flat under rescaling in the sense that \iyi1x^oo L{xy) / L{x) = 1 for every 
y > 0. For example, all powers and iterates of the logarithm are slowly 
varying at infinity. These are the class of admissible corrections to the 
power law x^^ . 

Part (i) of the theorem is an assertion of rigidity of scaling limits. We 
assume only that X{t) is measurable, positive and limt_»oo A(t) = oo. It then 
follows that the limits must define self-similar solutions, and A(t) must be 
the time scale associated to the self-similar solution, up to a slowly varying 
correction. Part (ii) of the theorem, and the sufficient conditions of Gallay 
and Mielke, show that the domains of attraction are determined by the tails 
of the initial data. As a consequence of part (i) of the theorem, the condition 
(1.11) is optimal. 

Gallay and Mielke [8] used a rather delicate Fourier analysis to establish 
the existence for self-similar solutions by studying their densities. We will 
use the proof of part (ii) of the theorem above to simplify much of this anal- 
ysis and extend it to the case when pfc 7^ for infinitely many k. In this we 
are motivated by a certain resemblance of the min-driven model to hydrody- 
namic limits of what are called A-coalescents in probability theory [2] , which 
are clustering processes involving arbitrarily many multiple collisions. We 
find a curious fact, namely, when ^p^k log k = oo there is no self-similar 
solution with finite mass (first moment). Still, the theorem above correctly 
describes the domains of attraction. Solutions with finite mass approach 
the self-similar solution with = 1, but this self-similar solution has infinite 
mass. 

1.4 Eternal solutions 

Theorem 1.1 is a particular example of the principle that the asymptotic 
behavior under rescaling is determined by the tail of the initial distribution. 
The dynamics exhibit sensitive dependence on initial conditions, as arbitrar- 
ily small changes in the tail of the initial data can lead to widely divergent 
asymptotic behavior. This indicates a kind of chaos, and it is of interest to 
find a precise formulation. A comprehensive analysis of such phenomena for 
the solvable cases of Smoluchowski's coagulation equation appears in [11]. 
This analysis is guided by an analogy with the probabilistic notion of infinite 
divisibility. Let us first describe these results informally. 

Clustering is an irreversible process, and in general we do not expect to be 
able to solve (1.3) backwards in time (for t < to). However, the self-similar 
solutions have the remarkable feature that they are defined for all t > 0. 
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That is, they are divisible under the coalescence process. We call a solution 
eternal if it is defined on the maximal interval (0, oo) consistent with (1.6). 
In probability theory, the infinitely divisible distributions are characterized 
by the celebrated Levy-Khintchine formula. In [11] we extended a result 
of Bertoin [1] showing that the class of eternal solutions to Smoluchowski's 
coagulation equations is also characterized by a Levy-Khintchine formula. 
Heuristically, this formula describes the emergence of eternal solutions from 
infinitesimally small clusters at i = 0. We also showed that the set of all 
subsequential limits — the scaling attractor — is in a one to one correspon- 
dence with the eternal solutions. A rigorous description of chaos is based 
on the fact that nonlinear dynamics on the scaling attractor is reduced to 
linear scaling dynamics using the Levy-Khintchine formula. 

In this article, we take the first step towards establishing a similar picture 
for min-driven clustering. Namely, we prove a Levy-Khintchine formula 
characterizing all eternal solutions for min driven clustering (Theorem 7.1). 
The choice of time scale (1.6) is very convenient for this analysis. 

1.5 Outline 

The rest of the article is organized as follows. We describe the solution 
procedure of Gallay and Mielke in the next section and discuss how the 
number-driven time scale is motivated by the important example of initial 
data that are monodisperse (a Dirac delta). The treatment here is formal. 
We establish some analytic prerequisites in Section 3. This is followed by 
rigorous results: the proof of well-posedness for measure- valued solutions 
(see Theorem 4.3) in Section 4, the study of self-similar profiles in Section 5, 
the characterization of domains of attraction in Section 6, and the charac- 
terization of eternal solutions in Section 7. 

2 The solution formula for min-driven clustering 

2.1 The generating function and moment identities 

As in the theory of branching processes, it is convenient to keep track of the 
clustering process with a generating function 

oo 

Q{z) = Y,Pkz'. (2.1) 

k=l 

For example, binary clustering as in the Allen-Cahn model corresponds to 
Q{z) = z^ . The generating function Q is analytic in the unit disk \\z\ < 1}, 
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and absolutely monotone (that is, Q and all its derivatives are positive on 
[0, 1)). We assume that the expected number of clusters in the mergers, 
denoted Qi, is finite. That is. 



Qi = Q'{1) = Y,kpk< 



oo. 



(2.2) 



fc=i 



We define a convolution operator 
rewrite (1.5) in the form 



^"kLiPkPt associated to Q, and 



dtfit,x) = f{t,l)i{Q{pt)ix-l)-Q^Pt{x)), x>/(t). (2.3) 

We extend the evolution equation (2.3) from densities to measures as 
follows. We consider a number measure I't and a probability measure Ft 
that are related to the densities (when they exist) by 



i^tidx) = f{t,x)dx, Ft{dx) 

Let M+ denote the interval [0, oo). If a : M+ 
support, then formally 

d 



ut{dx) 



pt{x)dx. 



(2.4) 



N{t) 

C is continuous with compact 



— / a{x)f{t,x)dx = / a{x)dtf{t,x)dx - a{l)f{t,l)l. 
We substitute for dtf{t,x) using (2.3), (2.4) to obtain the moment identity 



— / a{x)vt{dx) 



fc>i 



(2.5) 



/■ k \ k 

I' + ^yij - «(0 - X] '^^y^ 

\ i=l J i=l 



llFt{dy^). 



Some basic properties of the model are obtained by choosing suitable test 
functions a in (2.5). We set a{x) = x to see that mass is conserved: 



d_ 
di 



xut{dx) = 0. 



(2.6) 



When a = 1, we obtain the rate of change of the total number of clusters, 
N = -Qi/(t, /)/ = -NQipt{l)l. (2.7) 



1 + Y.y^] ""(^) 



llFtidyi). 



i=l 
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We substitute (2.7) in (2.3) to see that pt satisfies 

dtPt = ptil)iQipt)ix-l), x>l. (2.8) 

Similarly, we use (2.2), (2.5) and (2.7) to obtain the moment identity 

— / a{x)Ft{dx) = pt{l)l'S2pk I 

(2.9) 
This identity becomes an appropriate weak form of (1.7) after we later im- 
pose our choice of time scale t = 1/N, which means from (2.7) that 

Pt{l)i = 7^. (2.10) 

2.2 Gallay and Mielke's solution formula 

A remarkable feature of these min-driven clustering models is that the evolu- 
tion equation admits an elegant solution via the Fourier (or Laplace) trans- 
form. Our analysis relies heavily on this solution procedure, due to Gallay 
and Mielke [8]. The main difference with [8] is that we prefer to use the 
Laplace transform, denoted by 

Ptiq)= I e"^''ptix)dx, q>0. (2.11) 

Jr+ 

We set a{x) = e~'^^ in (2.9) to obtain the ordinary differential equation 

dtpt{q) = -iPtm e-'^' (1 - Qiptiq))) ■ (2.12) 

In order to integrate this equation, we define an analytic function (p via 

(This definition of <p differs by the factor Qi from that used in [8].) if is 
strictly increasing on [0, 1). In the case of binary clustering, Q{z) = z"^, and 
the functions (p and (p~^ are 

<^(z) = log (i±^) , ^-\w)=tanh^. (2.14) 

We substitute (2.13) in (2.12) to obtain 

dMpt{q)) = -{QiPtm e-"', q > 0. (2.15) 
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The choice of time scale has played no role in the analysis this far. Gallay 
and Mielke parametrize time by the minimum cluster size, setting 1 = 1. 
For clarity of notation, we denote this choice of time scale by r, reverting 
to the letter t when we introduce the number-driven time scale in (1.6). 

With 1{t) = T, the value of pr on the free boundary x = t plays an 
important role in the solution. We use this density to define a measure on 
(0, oo) that we call the trace measure A, with distribution function written 

.. N \ Ctn + [^ QiPs(s)ds, T > Tn, , , 

A{t) = { Jro^i^sV ^ ' - u, ^2.16) 

[qO T < To. 

Here tq denotes the initial time, and ag is any convenient constant. Equation 
(2.15) may now be rewritten 

rlA 

dMPM)) = -e-'^-r- (2-17) 

or 
Fix Ti > Tq. We integrate (2.17) from tq to ti to obtain 

viPrM) - viPrM) = - r e-'^'Aids). (2.18) 

J TO 

Since pr is supported in [r, oo), we have the estimate 

/oo 
Pr{x) dx = e-'^^ 

We now let ri — > oo in (2.18) to find the Laplace transform of A given by 

A{q)= [ e-^'A{ds) = ^{pr,). (2.19) 



Thus, the trace measure A is the inverse Laplace transform of fipro) a-i^d is 
determined completely by the initial data. 

We may now repeat this argument to determine the solution at any time 
T > Tq. We replace tq by r in (2.18), let ri —>■ oo, and obtain 

Pr{q) = ^'' {Mq)) , (2.20) 

where 

/•oo 

Ar{q) = / e-«M((is). (2.21) 
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We note that A^ is the Laplace transform of the truncated trace measure 
Ar satisfying 

Ar{ds) = H{s-T)A{ds), (2.22) 

where H is the Heaviside function. We may write 

Aris) = A{s) V A{t) (2.23) 

for the distribution function, with the notation a\/b = max(a, b). Therefore, 
the nonlinear evolution of pr is determined by the linear evolution of A^-. 
This global linearization underlies the analysis in [8]. 

2.3 The number-driven time scale and an extended solution 
formula 

It is natural to try to use the formula (2.20) as a basis for finding measure- 
valued solutions when the initial data pr^ is replaced by an arbitrary prob- 
ability distribution. However, we face two difficulties. The first is that it 
is not clear that (2.20) necessarily defines a measure pr- That is, it is not 
clear that the right-hand side of (2.20) is necessarily the Laplace transform 
of a measure. The second difficulty is that whenever the trace measure A 
has atoms, any solution defined through (2.20) is discontinuous in time, as 
is clear from (2.21). We overcome the first difficulty by an approximation 
argument. This relies on the simple and fundamental fact that a limit of 
completely monotone functions is completely monotone. We overcome the 
second difficulty by switching to the number-driven time scale (1.6). 

Henceforth, the letter t always denotes the number-driven time scale. 
The measure- valued solution is denoted by Ft, its Laplace transform by 



/>oo 

Ftiq)= / e-i^Ftidr) 
Jo 



and the minimum cluster size by l{t). These are related to solutions in the 
time scale r by 

l{t) = T, Ft{dx)=pr{x)dx, Ft{q) = pr{q). (2.24) 

We now rewrite the solution formula (2.20) in terms of measures. The 
relation l{t) = r, equation (2.7) and the definition of A in (2.16) imply 

— = Qip^{T)dT = A{dT). (2.25) 
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The differential equation (2.12) now takes the form 

-ql{t) 

dtFtiq) = -^5— (1 - QiFtiq)). (2.26) 

Wit 

If to denotes the initial time, we may integrate equation (2.25) to obtain 

log(^)=^(/(t))-ao. (2.27) 

The change of variables (2.24) and (2.25) affects the Laplace transform of 
Ar as follows: 

/■OO 7 

^-1^ A(^'r.\ = / ^-9'W"'^ 



A{dr)= / e-«'^^^ — . (2.28) 

Jt s 

These calculations yield the following revised solution procedure: Given 

an arbitrary initial probability measure F^g, the trace measure A is found 

as in (2.19) by inverting its Laplace transform, given by 

A{q) = ^{FtM)- (2-29) 

Next, we determine l{t) through inverting (2.27). Once l{t) is known, the 
solution Ft is determined by inverting the Laplace transform given as in 
(2.20) and (2.21) by 



Ft{<D = V' 



^i f r ^~<iiis)^\ t>to. (2.30) 



The main observation is that working with l{t) instead of A^ yields an 
evolution continuous in time. Since l{t) is an increasing function, it has at 
worst jump discontinuities. But (2.30) shows that discontinuities in / do not 
affect the continuity in t of Ft{q), and thus the continuity of Ft in the weak 
topology. 

2.4 An example: monodisperse initial data 

Let us illustrate the meaning of the extended solution formula in the new 
time scale with an important example. Set to = tq = 1, Q{z) = z'^ and 
Fi{x) = lx>i- That is, initially all clusters have size 1. Then Fi{q) = e~'^, 
and 

A{q) = V5(Fi(g)) = log(l + e-1) - log(l - e"-?). 

We differentiate with respect to q and simplify to obtain 

-d,A{q) = ^^^ = 2 (e-« + e'^^ + e~''^ + ....) 
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Since e '' is the Laplace transform of 6k{dx), the trace is 

A{x)= Y: I- (2-31) 

k<x, k odd 

A has jump discontinuities at the odd integers. We shall work with the right 
continuous inverse, so that the minimum cluster size is 

l{t)=k, te[tk_2,tk), tfc = e2(i+3+-i), A; odd, (2.32) 

with t-i = to = 1. The solution formula (2.30) now yields 

/ ti.\ e'^'i e"«(''+2) e-'?(^+'*) 

^(«(,)) = (^,„g-j _ + __ + __ + ..., ,,fc_„«. 

The solution has the following interpretation. For t G \tk-2itk)-, Ft is sup- 
ported on the odd integers greater than or equal to k. The fraction of clusters 
of size k decays continuously to zero over the time interval \tk-2-,tk)- Thus, 
the number-driven time scale regularizes jumps in the trace measure by pro- 
viding a finite time for these jumps to vanish. A moment's reflection suggests 
that this is what we should expect if we approximate the monodisperse data 
by a smooth density. 

3 Analytic preliminaries 

This section is a summary of the main analytic methods we use. We present 
some facts about distribution functions, Laplace transforms and Tauberian 
theorems in a form suitable for use in later sections. 

3.1 Distribution functions 

We will consider measures on an interval J C M. We study a measure 
through its distribution function, often using the same notation for both. It 
is therefore convenient to introduce the following conventions for brevity. A 
distribution function /: J ^ M is a right-continuous, increasing function. 
("Increasing" means that xi < X2 implies f{xi) < f{x2)-) A distribution 
function / is identified with a measure via 

/((X1,X2]) = /(X2)-/(X1). (3.1) 

We do not assume that the function / is positive, since the trace measure 
for self-similar solutions is of the form A{t) = 01ogr, r > 0. Following 
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probabilistic convention, we say a sequence /„ of measures on J converges 
weakly to / (written /„ -^ f) if and only if /n([a, b]) -^ f{[a,b]) as n ^ oo 
whenever a,b £ J are not atoms of /, meaning /({a}) = /({&}) = 0. We 
have fn^f'^^ and only if at every point of continuity of f{x), 

fn{x) + c„ ^ /(x) as n ^ oo, 

for some constants c„ independent of x. 

Given a distribution function / : J ^ M, its epigraph is the set 

r(/) = {{x,y) e R2 I f{x-) <y< f{x),xe J}. 

There is a unique distribution function f', with epigraph obtained by reflec- 
tion through x = y: r(/1') = {{x,y) \ {y,x) G r(/)}. We call /^ the inverse 
of /. We can write 

/t(r)=inf{tG J|/(t)>r}, r < sup{/(t) | t G J}. (3.2) 

The following convergence result is not difficult to prove. 

Lemma 3.1. Suppose fn'- J ^ M is a sequence of distribution functions that 
converges to f at all points of continuity. Then for every point of continuity 
X of f\ fn{x) is defined for sufficiently large n, and lim.n^oo fJt{x) = f^x). 

Lemma 3.2. Suppose f: J ^ M is a distribution function. There exist 
monotonically increasing and monotonically decreasing sequences of piece- 
wise constant distribution functions that converge to f at all points of con- 
tinuity. 

Proof. We assume that J = [0, co) for clarity, and we only construct an 
increasing sequence of approximations. The argument is easily generalized. 
Suppose /c > 1 is an integer. Let Ifc denote the indicator function lk{x) = 1, 
fc < X < A; + 1, and Ifc(x) = otherwise. Consider the sequence of functions 

9n{x) = X]-^ ( ~) '^kinx). 

Since / is increasing, we have 

gi{x) < g2{x) < ... < gn{x) < f{x). 

Let us establish convergence of gn- Given x > and n, let /c„ denote the 
integer such that kn < nx < kn + I. Clearly lim„_»oo hi/n = x. Suppose 
X is a point of continuity of /. Then gn{x) = f{kn/n), therefore gn{x) -^ 
fix). D 
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3.2 Min history and trace 

Our analysis will focus on two related distribution functions — the minimum 
cluster-size history / and the trace A. 

Fix to > 0. A min history is a positive distribution function / on [fg, oo). 
In this article we also require a min history to be unbounded: /(t) ^ oo as 
t ^ oo. Given any min history /, we associate a trace A on (0, oo) via 

^(r) = log/^(r), r>0. (3.3) 

Note A{t) = log to for < r < /(to) due to (3.2). 

Conversely, we say ^ is a trace if it is a distribution function on (0, oo) 
such that (i) A{t) = log to on some nonempty, maximal interval (0, tq), and 
(ii) A{t) ^ do as r ^ oo. Given any trace A, we can associate a min history 
^by 

l{t) = expA^{t), t>to. (3.4) 

Proposition 3.3 (Change of variables), (a) Assume I is a min history and 
there exists c > such that 

oo Jj. 

e-«'W— <oo, gG(c,oo). (3.5) 

Then the trace A given by (3.3) satisfies 

l-CO f-CO Jj. 

A{q)= e-'^^A{dT)= e~«^W — , gG(c,oo). (3.6) 

Jo Jto * 

(b) Assume A is a trace such that A{q) < oo for q G (c, oo), and let l{t) be 
defined by (3.4). Then (3.6) holds. 

Proof. 1. We first verify the equality for piecewise constant functions. Sup- 
pose < to < ti < t2 < . . . and Q < Iq < li < . . . are strictly increasing 
sequences. Consider the piecewise constant, increasing function 

oo 

m = Y,h\t,,t,^,){t), t>to. (3.7) 

fc=0 

An associated trace is given by 

oo 

A{t) = Y, logtk+ilii,,i,^,){T), r > 0. (3.8) 

k=-l 
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Here /_i = 0. Conversely, if A is defined by (3.8), then / is given by (3.7). 
We fix (7 G (c, oo), assume (3.5), and use (3.7) and (3.8) to compute the 
integrals in (3.6). Both equal 

-g'o ,^„ /'M , ,-g«- /^*2\ , _„,,, ft3 



e-9'o log / ^ 1 + g-g'i log / ^ 1 + e-9«2 log / ^ l + . . . 

2. Suppose / is given, and (3.5) holds. Fix g > 0. We approximate 
/ by a decreasing sequence of piecewise constant functions Z„ J, L Then 
Q-qin I g-qi gince q > c > 0, and moreover we have An j ^ by Lemma 3.1. 
By the monotone convergence theorem, 



oo 
3 


e" 


-ql{t) 


dt 
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= hm 

n^ooj 
■oo 


/•oo 
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qlr. 
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■oo 
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e-'^^A 


-nidr) 
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/ 





r«M(dr). 

"-"^JO JO 

3. To prove (b), assume A is such that A{q) < oo for (^ € (c, oo). Consider 
a sequence of piecewise constant increasing functions An ] A. Then /„ J, I 
and we may apply the monotone convergence theorem again. D 

Proposition 3.3 and Lemma 3.1 allow us to reformulate the classical 
equivalence between weak convergence of measures and pointwise conver- 
gence of Laplace transforms. The following theorem is a slight modification 
of [7, XIII.1.2a]. 

Theorem 3.4. Suppose In is a sequence of min histories on [tQ,oo) with 
associated traces An that satisfy 

snp An{q) < oo, g G (c, oo), (3.9) 

n 

for some c > 0. Then there is a min history I such that ln{t) -^ l{t) as 
n ^ oo at all points of continuity, if and only if there is a trace A associated 
to I such that An{q) -^ A{q) as n ^ oo, for all q G (c, oo). 

Proof. Suppose /„ -^ I. Lemma 3.1 and definition (3.3) then imply the distri- 
bution functions An -^ A. By the classical criterion for weak convergence of 
measures [7, Xin.l.2a], under the hypothesis (3.9) it follows An{q) —>■ A{q) 
for all q > c. 

Conversely, suppose An{q) -^ A{q) for all q > c, where A is a trace (in 
particular A{t) -^ oo as t ^ oo). By the classical criterion, the distribution 
functions An -^ A, and Lemma 3.1 yields In -^ I where / is given by (3.4). D 



18 G. Menon, B. Niethaninier and R. L. Pego 

3.3 Regular variation 

A measurable function R : (0, oo) -^ (0, oo) is regularly varying at infinity 
with index 6* G M (written R € RVe) if 

R(Xx) Or 

lim = X , for every x > 0. (3.10) 

In this case, R{x) = x L{x) where L is slowly varying at infinity. 

The class of regularly varying functions is remarkably rigid. For example, 
there is no need to assume that the limit in (3.10) exists for every x > 0. If i? 
is a positive, measurable function on the half-line and \mi\^oo R{\x) / R{\) 
exists, is positive and finite for x in a set of positive measure, then / is 
regularly varying at infinity with some index ^ G M. This fundamental 
rigidity lemma (see [7, VIII. 8.1] and [3, 1.4.1]) plays a key role in our analysis. 

The class RVg is of fundamental utility in Tauberian arguments linking 
a measure u on [0, oo) and its Laplace transform u{q) = L e~^^u{dx)^ see 
[7, XIII.5.2]: 

Theorem 3.5. If L is slowly varying at infinity and < 9 < oo, then the 
following are equivalent: 

i/(x) ~ X L{x), X ^ oo, (3-11) 

and 

iy{q)r^q-'L{l/q)r{l + e), q ^ 0. (3.12) 

Moreover, this equivalence remains true when we interchange the roles of 
the origin and infinity, namely when x ^ and q ^ oo. 

A refinement of this result will prove useful for us — de Haan's exponential 
Tauberian theorem [3, Thm. 3.9.3], see [5]. 

Theorem 3.6. exp z^(x) is regularly varying at infinity with index 9 if and 
only if ex];)i>{q) is regularly varying at zero with index 9. If either holds, 
then 

i^il/q) - v{q) ^j9, q^ 0, 
where 7 = 0.577215665 . . . is the Euler-Mascheroni constant. 
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3.4 Rate of divergence in the solution formula 

For several reasons, we need to study carefully the asymptotic behavior of 
(p{z) as z — > 1. For this, it is convenient to define K{q) via 

fl—q / f) 1 \ 

-^og.{q)=^{l-q)+logq = J^ (^^-^ - ^ d. (3.13) 

(note — log K{q) > 0) and set 

Ko := lim K(g) G [0, 1] . (3.14) 

q-^O 

With this notation kq corresponds to the number k introduced in [8]. For 
Q{z) = z^ we have kq = \- 

We will show that kq > if and only if ^^^=1^ log A;)pfc < co. Since in [8] 
the function Q is a polynomial this finiteness condition is always satisfied. 
However, we can characterize self-similar solutions and their domains of 
attraction also in the case kq = 0. 

Lemma 3.7. The junction K{q) as defined in (3.13) satisfies K{q) ^ as 
q —> Q if and only if Yl^=i Pkk log k = oo. 

Proof. We compute 

oo k—l 



B(-) := «. - if^ = E« (^ - ^) = Ek Ed 

k=l ^ ^ k=\ 3=1 

Note -R(l) = and the integrand in (3.13) is 

R{z) R{z) 1 R{z) 1 



z^). (3.15) 



l-Q{z) \-z Qi-R{z) 1-z Qi + o(l) 



(3.16) 



as z ^ 1 . Thus it suffices to show L R{z) dz/{l — z) < oo if and only if 
"^ Pkk log k < oo. We observe 



fe-i , i fc-i j , fe-ifc-i, fc-i 



iE^^-EEy-EEy-E^-^'°^* i^^") 



j=l ■ j=l 1=1 ' 1=1 j=l ' 1=1 



as k ^ oo, whence the desired result follows from (3.15). D 

Lemma 3.8. The function K{q) is slowly varying as q ^> 0. 
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Proof. We need to show that — log{K{Xq) / K{q)) ^ as g ^ for any A > 0. 
From (3.15) and (3.16) above, we see that if we replace the integrand in 
(3.13) by R{z)/{1 — z), then in integrating from 1 — Xq to 1 — q the error is 
only o(l) as g' ^ 0. Hence by (3.15), 

-.o«^^i:..^E(^^-V^)-a. 

^^^ k=i 1=1 
as g' ^ 0. We consider without loss of generality A > 1 and estimate 

On the other hand 

^ I I 

1=1 

as q ^ for any k and thus the claim follows from the dominated conver- 
gence theorem. D 

4 Well posedness for measures 

In this section, we first work with the min-driven time scale 1{t) = r used 
by Gallay and Mielke. We prove that a continuous probability density with 
min To > defines a solution to (2.8). This is a weaker form of the well- 
posedness theorem of [8]. It is included for its simplicity, and because it 
is the basis for weak solutions. We then switch to the number-driven time 
scale (1.6) and use the moment identity (2.9) to show that (2.8) defines a 
continuous dynamical system on the space of probability measures P(M+). 

4.1 Classical solutions 

The solution formulas of section 2.2, while explicit, are not immediately 
suited for the construction of solutions. The main difficulty is to show that 
positive initial data yields a positive solution. We construct solutions by 
rewriting (2.8) in integral form 



pr{x) = p{x)+ ps{s)Q{ps)ix - s)ds, X > T, (4.1) 

J TO 

Pt{x) = 0, X <T. 

We then have 
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Theorem 4.1. Suppose p is a continuous probability density with positive 
min tq. There exists a unique solution to (4-1) on [tq, oo) such that p^-o = P 
and the solution has the following properties. 

(a) For every t > tq, pr is a continuous probability density with min r. 

(b) The solution formula (2.20) holds for every (7 > 0, t >tq. 

Proof. We sketch a proof of existence similar to the direct approach in [4] 
for Q{z) = z^ using a different time scale. We fix tq > and let p be given. 
Note that since the solution is to satisfy Pt{x) = for x < r, the convolution 
term on the right-hand side of (4.1) will depend only upon values of Priv) 
for Tq < y < X — T < X — Tq. In particular, this convolution term vanishes 
for X < 2tq. 

This means we can construct the solution for tq < r < 2tq by an in- 
ductive procedure as follows: For tq < t < x < 2ro we have Prix) = p{x) 
and in particular Prij) = p{t). For tq < r < 2ro, successively on strips 
X £ [fcro, {k + 2)ro), for k = 2,4, . . ., by simple integration in time we can 
now compute Pt{x) from (4.1), where the right-hand side is always known 
from a previous step. This determines Pr{x) for tq < r < 2tq and all x. 

To determine the solution globally for all r > tq, the idea is to replace 
2tq by To and repeat. But in order to justify this we need to verify that p-r 
remains integrable and conserves total probability. In particular we need to 
justify (2.12). Let us introduce the distribution function 

t-x nx 

Rr{x)= p^{y)dy= Priy)dy. (4.2) 

This is the probability that a domain has size < x at time t. Note that for 
any two distribution functions R{x), R{x) on [0, 00) we have 

px px 

R*R{x)= / R{x-y)R{dy)< j R{x)R{dy) = R{x)R{x). 
Jo Jo 

Integrating the convolution term in the integrand of (4.1), we find 



00 CO 



J Q{Pr){y -r)dy = Y^PkKH^ -r)< Y,PkRr{xf = Q{Rr{x)). 



fe=l k=l 



Then it follows drRrix) < pr{T){Q{Rr{x)) — 1), and since pr(j) > and 
Rrix) < 1 initially, Rt-{x) is decreasing in r for fixed x. It follows Rr {00) < 
1, and so the Laplace transform 

Rr{q)= e-i''Rridx)= e-«^/)^(x)dx 
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is well defined and < e^''^. Since drRri-c) is continuous in r for all x, Rriq) 
is C^ in r for all q > 0. This justifies (2.12) and the the computations leading 
up to (2.19) and the solution formula (2.20). From (2.19) we deduce that 
since i?ro(0) = 1, ^(0+) = oo = 1^(0+) and then (2.20) yields Rr{0^) = 
1 = Rr{oo), proving that pr is a probability density for all r. D 

The next lemma provides a uniform estimate for smooth approximations. 

Lemma 4.2. Suppose p is a continuous probability density with positive niin 
tq and pj- is the solution to (4-1) with initial data p. Then (with a^ = in 
(2.16)), 

^(r)<y%logf-). (4.3) 

log 2 V ^0 / 



Consequently, 



l{t) 1 /t\(l°s2)/Qi 



/(to) " 2 \tQ 



Proof. If m is an integer such that r € [2™^-'^ro, 2™'ro) we divide the domain 
of integration [ro,r) into m pieces and use the fact that ps[s) = Pr{s) for 
To V gS < r < s to obtain 

Q-^A{t) = / prf,{s)ds+ P2To(s)ds + ... / p2m-i^^^{s) ds 

Jto J2to J2™-1to 

^ ^ log (r /to) 

< m < — ; h 1. 

- log2 



D 



4.2 Weak solutions 



We now switch to the number-driven time scale N{t) = 1/t. In order to de- 
fine weak solutions, we fix a test function a, substitute (2.10) in the moment 
identity (2.9) and integrate in time between to and t to obtain 



a{x)Ft{dx) - 


1 a{x)Ft^,{dx) -- 


„t oo „ 


r / k 



k 


{dVi) 


ds 
Qis' 



a{l{.s)) 



We will consider continuous test functions with lim^;— >oo o,{x) = 0. Let 
Co(M_|_) denote the space of such functions with the topology of uniform 
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convergence. Let "/^(M-i-) denote the space of probability measures on M_|_ 
equipped with the weak topology. Assume to > is fixed. 

Definition 4.1. Let J C (0, oo) be an interval. We say that a map F : J ^ 
■p(R+) is a weak solution for min-driven clustering on J if 

1. The map i '-^ Jj^ a{x)Ft{dx) is measurable for every a € Co(M_|_). 

2. The min of F(, denoted l{t), is positive and increasing. 

3. The moment identity (4.5) holds for each a G Co(M4,) and t,to £ J- 

Theorem 4.3. (a) Suppose F G V{^j^) has positive min, and to > 0. 
Then there is a weak solution F for min-driven clustering on [to,oo) 
with FfQ = F. Moreover, the min l{t) of Ft satisfies (4-4) ■ 

(h) The solution in (a) is unique on [to,ti] for any ti > t^. 

(c) Let F^"^' be a sequence in Vi^^) with positive min and F^^' the weak 
solutions with F^ = F'"\ Assume liran^oo F^^) = F and the limit 
has positive min. Then F^ -^ Ft for every t > t^. 

Proof. 1. It follows from Weierstrass' approximation theorem that finite 
linear combinations X]fc=i c^e" are dense in Co(M+). Therefore, in order 
to verify the moment identity, it is sufficient to consider the test functions 
Q-qx ^ g > 0. Thus, to prove existence of a weak solution on [to, oo), it suffices 
to construct F weakly continuous such that the Laplace transform satisfies 
the solution formula (2.30). 

2. Let To denote the min of F. We approximate F by a sequence of 
continuous probability densities p'""' with min Tq with lim„^ooTo = tq. 
We further assume that p^"^' is strictly positive on [tq , cxd). It is immediate 
from (4.1) that the solutions pf" are strictly positive on [r, cxd). The trace 
for these solutions, A^ "i is obtained from (2.16) with tq replaced by Tq , p by 
p'™^ and the choice uq = log to. A^"^' is continuous and strictly increasing, 
thus so are the min histories /'-"'•*. We change variables from the solution 
formula (2.20) to (2.30) to obtain 



^-(n)(,) = i^W(,) = ^~ifre-'^''"'W^V t>to. (4.6) 



3. As n ^ oo, we have p^"''{q) -^ F{q) for every q > and < F{q) < 
1 for q > 0. The measures A^""' are supported on [tq ,oo) and satisfy 
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^(«)(g) = ^{F}^\q)). Therefore, 

lim A'-''\q)=ip{F{q)), q > 0. 

It follows that the traces ^'"' converge weakly to a trace A supported on 
[ro,oo) that satisfies A[q) = (p[F{q)). Therefore, by Theorem 3.4 the min 
histories /*•"-* converge to /, the min history associated to A. Since A"' satisfy 
the uniform estimate (4.4), we may use the dominated convergence theorem 
to assert 

limF(")(,) = V.-ifre-^'W^V t>t,. 

This shows that for every t > to the measures F^' converge weakly to a 
measure Ft that satisfies (2.30). Then Ft is a probability measure since 
Ft(0+) = f^^Hoo) = 1- It similarly follows from (2.30) that Ft -^ Ft^ as 
t —> ti for every ti G [to, oo). This completes the proof of part (a), except 
that it remains to show that l{t) is in fact the min of Ft. 

4- For (b) it suffices to prove uniqueness on [tojti] for some ti > to- The 
key is to prove uniqueness of the min l(t) on such a time interval, since then 
uniqueness of Ft is easy to establish via the Laplace transform. Note that 
any weak solution is weakly continuous in time, since the right-hand side of 
(4.5) is Lipschitz continuous for any test function a G C^(W-^-). Then since 
F(q(x) > for all X > tq, the min l{t) is right continuous at to, so there 
exists ti > to with tq < ti = l{ti) < 2tq. Now the idea is that for clusters 
of size less than 2ro, there is no gain, only loss. We claim that 

Ft{x) = (Ft,{x) -l-iog(^\\ , t G [to,ti], X G [ro,2ro). (4.7) 

This follows by first considering points x of continuity of both Ft and Ft^ 
such that l{t) < X < 2ro, so Ft{x) > — approximate l[o,a;) by continuous 
test functions a supported in [0, 2ro). By right continuity and the definition 
of min, (4.7) follows for ah x G [/(t),2ro), and Ft{x) = for x < l{t). Now 
by weak continuity in time, we infer that with 

i{x) = toexp((5iF(x)), 

for each point of continuity of F in (tq, ti] we have l{t) < x for t < t{x), and 
X < l{t) for t{x) < t. Now clearly l(t) is uniquely determined by F, since it 
is locally the inverse of t. It follows from taking the Laplace transform that 
the min history l{t) constructed in (a), that satisfies (2.30), agrees with the 
min of Ft- 

5. Part (c) is proven by an argument very similar to Step 2 above. D 
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Scaling. We note for use below the following scaling property that follows 
easily from the moment identity (4.5): Let a,b > 0. If F is a weak solution 
for min-driven clustering on an interval J, then F is a weak solution on J/a, 
where 

Ft{x) = Fat{bx), teJ/a,x>0. (4.8) 

5 Self-similar solutions 

Let us recall that Gallay and Mielke have classified the self-similar solutions 
as in (1.9) in the case when pk = for all large k. Here we will recover and 
extend the basic existence results and classify all domains of attraction by 
simple means based on the Laplace transform and Tauberian arguments. 

Existence. In terms of the distribution function Ft, self-similar solutions 
are of the form Ft{x) = F^:{x/l{t)) for some distribution function F^ with 
positive min (= 1) and some min history l{t). Without loss of generality we 
can assume t^ = 1 and /(to) = 1- From (2.30) we obtain for the Laplace 
transform of F^ that 

(5.1) 
Hence we conclude that the min history of self-similar solutions must satisfy 
l{ts)/l{t) = g{s) for some function g{s). Since I is also positive, increasing 
and non-constant, we conclude (see [7, VIII.8.1], or subsection 6.1 below) 
that necessarily, for some ^ > 0, l{s) = s^'^ and that F* must be a distribu- 
tion F^ ' that satisfies 

F('Hq) = ^''ieEi{q)), (5.2) 

where Ei(g) = J e^^ds/s denotes the exponential integral. Provided such 

a distribution F(^) does exist, then Ft{q) = F^^\t^/^q) satisfies (2.30) and 
is continuous in time, hence determines a self-similar solution by step 1 of 
the proof of Theorem 4.3. The corresponding trace measures are given by 

A(^) (r) = 6* log r, r G (0, oo) . (5.3) 

As a byproduct of the characterization of scaling limits in section 6, we will 
prove the existence of distributions F^^' that satisfy (5.2) for < ^ < 1, and 
that 9 < 1 is necessary. 
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Densities. Next we show that for the self-similar solutions F^ \x/t^' ) 
{0 < < 1), the probability distributions F^^' have piecewise smooth den- 
sities p^ ' that satisfy an integrodifferential equation, 

-a,(ypW(y)) = ^Q(pW)(y-l), y > 1, pW(l) = i-, (5.4) 

and we study the decay of p^ '{y) as y — > 00. When Pfc = for all large k (so 
Q{z) is a polynomial) this has been done by Gallay and Mielke [8]. We will 
recover most of their results by simpler means (the exponential rate of decay 
for = 1 is an exception) and extend them to the case when pk is nonzero 
for infinitely many k. It turns out, however, that when Y^p^klogk = 00, 
none of the profiles have finite first moment, including the case 9 = 1. 

Let 6 G (0,1]. We claim that the probability distribution F = F'^^^ 
satisfies the following weak-form profile equation: 



fe \ \ k 



(5.5) 
for all (7^ functions a € Co(M-i-). This follows from the fact that we know 
F[x/t^' ) is a weak solution for the min-driven clustering equation (4.5) 
with l{t) = t^' . Changing variables in (4.5) via x = l{t)x and similarly for 
j/j, we differentiate at t = 1 to obtain (5.5). 

Now, the min of F is 1 (this will be shown in section 6.2), so F has 
density p{x) = on (0, 1). Taking a to be supported in (1, 2) we find that 
the right-hand side of (5.5) vanishes. Hence restricted to (1, 2), the measure 
xF{dx) = (3 dx for some constant /3, so F has density p{x) = (5/x on (1, 2). 
Taking a{x) = for x < 1, a{x) = 1 for x > 2 (approximated by limits — 
note J" xF{dx) -^ along a subsequence), we find [3 = 6/Qi. Taking 
a supported in (0,2) , we find the right-hand side is — /3a(l) and we can 
conclude that 1 is not an atom of F. 

In this way, proceeding inductively on intervals (0, n) we deduce that 
F has density p{x) satisfying (5.4). When a has support in (l,n + 1) the 
right-hand side depends on the restriction of F to (0, n) where it has density 
p{x)., and thus xF{dx) has density xp{x) determined by (5.4) on (l,n + 1). 

Decay. Next we wish to characterize the decay behavior of the densities 
p(^). For this we use the properties of the function K{q) which was introduced 
in (3.13). We denote by k* the de Bruijn conjugate of k. See [3, Sec. 1.5.7]. 
This is a slowly varying function satisfying 

K{q)K {qK{q)) ~ 1 as q ^ 0. (5.6) 
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If kq > 0, then k*(0^) = Kq ^, and if kq = 0, then k*(0'*') = oo. 

Theorem 5.1. For every 9 G (0,1], the density p^^'{x) of the self-similar 
profile F^ ' has the following properties: 



(i) If £ (0, 1), then as x —> oo, 



e{i-e) 



p(^)(x)~x-(i+^)e^-K#(.-^)^^. (5.7) 

Here 7 is the Euler-Mascheroni constant and T is the T-function. 
(a) If 9 = 1 then as x ^ 00, 

yp(^)(y)dy ~e^K#(x"^) (5.8) 



Remarks: The asymptotics (5.8) imply the result of [8] for total mass in 
the case that kq > 0. We will not pursue here the delicate question of the 
precise (exponential) decay rate of the density p^^' in this case, which was 
studied in [4] for Q{z) = z^ and in the polynomial case in [8]. If kq = 0, 
however, we see that F^^> has infinite mass. 

Proof. We rewrite equation (5.2) for F^^> using (3.13) and standard asymp- 
totics for the exponential integral [14], as follows. With w = 1 — F^ >{q), 

-^{F^^\q)) = \og{wn{w)) = -0Ei(g) = ^(logg + 7 + o(l)), g ^ 0. 

Then wk{w) ~ q e '^ , whence asymptotic inversion [3, Sec. 1.5.7] yields 

w = l-F^^\q)r^q^e^^K*{q^). (5.9) 

Now differentiating (i.e., using Lemma 3.3 of [10]) we find 

POO 

/ e-«^xF(^)(dx) = -dgF^^Hq) ~ eq^-^e^''K,*{q^), q ^ 0. 

Jo 

Now the Tauberian theorem 3.5 implies 

yF^'Hdy) ~ ^J x^-'e'^K*{x'') , x ^ 00, (5.10) 







which is just (5.8) in the case 6 = 1. 

We can now also derive the decay behavior of p^ > in the case 9 < I. In 
fact, it follows from (5.4) that yp '{y) is decreasing. Hence by a lemma in 
[7, XIII.5], (5.10) implies (5.7). D 
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6 Domains of attraction of self-similar solutions 

We proceed to prove Theorem 1.1. The proof has two parts. The first is to 
show that regular variation of l{t) as t ^ cx) is necessary and sufficient for 
convergence to a scaling limit. The second is the equivalence between regular 
variation of l{t) as t ^ oo and regular variation of J^ yFtf^{dy) as x ^ oo. 
The second part is based on the Tauberian theorems 3.5 and 3.6. The most 
subtle aspect (despite the simple proof) is to deduce regular variation of / 
from the existence of a scaling limit. This is the assertion of rigidity, and 
we treat it first. 

6.1 Regular variation of the min history is necessary 

1. Assume there is a rescaling A(t) -^ oo and a probability distribution 
function F* with positive min (called r*) such that Ft{\{t)-) -^ F^. This is 
equivalent to convergence of the Laplace transforms, 

After a trivial scaling of time and cluster size, we may assume to = 1 and 
r* = 1. Since F^, is a probability measure with positive min, there is a unique 
trace measure A^:, with A^(t) = on (0, 1) and A*(r) ^ oo as r ^ oo, and 
a min history l^ on [1, oo), with ^^,(1) = 1, such that 

AM= / e~'^^*(^)- = ^{F4q)). (6.2) 

Now, s 1-^ FtsiX{t)x) is a rescaled solution, by (4.8). The rescaled min 
histories given by l^^>{s) = l{ts)/X{t), t,s > I have associated trace measures 
A^^' satisfying A^^'{q) — > A^,{q) as t ^ oo, for all q > 0. We use the solution 
formula (2.30), (6.1) and (6.2) to obtain 

i^ r^-'^m/xit)^^ r^-,i4s)^^ ^>o. (6.3) 

t^coj^ S Ji S 

Theorem 3.4 now implies that at every point of continuity of /*, 

lim l{ts)/X{t) = h{s). (6.4) 

t— >oo 

2. Let So be a point of continuity of Z*, and let U{t) = l{tsQ), t > 1. 
Let B be the set of x > 1 such that 

^|;{x) = hm ^ (6.5) 

t-*oo U (t) 
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exists. We deduce ^(x) is a power of x by following the simple argument in 
[7, VIII.8.1]: By (6.4), we have x G i? if xsq is a point of continuity of l^. If 
xi, X2 £ B then xiX2 € B and 

i^{xiX2) = ^l^{xl)^p{x2). (6.6) 

Since U is increasing, so is ip, and since B is dense in [1, oo) we can extend 
ip by right continuity so (6.6) holds for all xi,X2 > 1. Then we can set 
ip{x) = l/^(l/x) for X G (0, 1) and have (6.6) for all x > 0. Since ^(x) is 
positive, locally bounded and not constant, it follows ^ is a pure power law, 
and we can write tp{x) = x^' for some ^ > 0. 

3. Since ip is continuous and increasing, it is easy to see (6.5) holds 
for all X G [l,oo) = B. Then we infer L(t) = U{t)t^^' is slowly varying, 
and U , hence I, is regularly varying at cx) with index 1/0. Further, since 
1*{xsq)/U{sq) = il){x) = x^'^ whenever xsq is a point of continuity of /*, and 
So is an arbitrary point of continuity, it follows h{t) = t^'^ for alH > 1. 

6.2 Regular variation of the min history is sufficient 

Let us assume that the min history / is regularly varying, as in (1.12). 
Convergence to self-similar form is then quick: Since 



s 



m L{t) 

we have 

hm iM = ,1/^ 

Moreover, for any e > 0, there exists t^ > such that 

i/e-e < iM t > t 

' - l{t) ' -^'- 

Therefore, we may use the solution formula and the dominated convergence 
theorem to see that as t ^ oo 

^ (Pt (^)) = re-'(*^)/'(*)- - fe--^' -. (6.7) 



Kt)JJ J I s Ji s 

This implies limj^oo -^t(K^)^) = F^'{x) for every x > 0, where F^' is 
a probability distribution with positive min (= 1) and Laplace transform 
given by (5.2). 
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6.3 Tauberian arguments 

We now prove that if < < 1 and J„ yFt^^{dy) is regularly varying with 
index 1 — as x -^ oo, then l(t) is regularly varying with index 1/6 as t ^ oo, 
and conversely. 

1. It is convenient first to assume (1-11) in the form 

a 

yFt^,{dy) - — - — — x^~^Li(x), x ^ oo, 
il^-^j 

where Li is slowly varying. Then Theorem 3.5 implies 

-dgFt,{q)^9q'-'Li{q-^), q ^ 0. 
By integration (not difficult to justify as in Lemma 3.3 of [10]) we find 
w:=l-FtM-Q'Li{q~^), q ^ 0. 

Using the solution formula (2.30) and (3.13), we write 

/■oo ^ 

log {wkH) = M^toiq)) = - / e-'?'(^)- = -A{q) . 

Jt s 

Hence, exp{—A{q)) is regularly varying with exponent 9, and Theorem 3.6 
(de Haan's exponential Tauberian theorem) implies that 

^Ail/q)^-A{q) ^ /t(g-l)y;^(y;) ^ g^^ ^ g ^ Q. (6.8) 

Thus, /t(T) ~ eT^r^Li(T^) as r ^ oo, with 

1 



LUs) 



Li(sVe)K(s-iLi(si/e))- 



It is easy to show Li is slowly varying, using the uniform convergence theo- 
rem for slowly varying functions [3, Theorem 1.2.1]. Asymptotically solving 
t = 1'{t) by inverting s i-^ sLi{s) using the de Bruijn conjugate [3, Theorem 
1.5.13] finally yields 

l{t) ~ e-^t^/^L* (t)l/^ t ^ oo , (6.9) 

which gives (1.12). If kq := liniq^o f^iQ) > 0, we also can write 

l{t) ~ e-^'{Kot)^/\L^'^^'')*{t), t -^ oo. (6.10) 
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2. We now prove the converse. Assuming that (1.12) holds, then 

l{t) ~ e-^tl/^L2(t)^/^ t -^ oo, (6.11) 

for some slowly varying function L2. Then, by inversion, /^(r) ~ e^^ t^ L2 {t^) 
as r — > 00. Since 1\t) = expj4(T) and A{q) = —\og{wK{w)) with w = 
1 — Ft^{q), we infer from Theorem 3.6 that (6.8) is true. This implies 

wk{w) ^q^L*{q-'^)-\ q ^ 0, 

and asymptotic inversion yields w ~ q^L2{q^^) as g ^ 0, with 

and L2 is slowly varying. Differentiating using Lemma 3.3 of [10], we find 

dgW = -d,Ft,{q) ~ 6q'-'L2 (g"^) , q ^ 0. 
Since —dgFtg (q) = L e^'^^xFt^ (dx) the Tauberian Theorem 3.5 implies that 

^ yFtMy)-^ ^^2^_^-^ ^'^'L2{x), x^cx), (6.12) 

which is just (1.11). If kq := lim^^o f^iQ) > 0, then 

/ vFtn(dy) r-^ —-, T^r^-iT , X — > 00 . (6.13) 

J^ tov yj n2-e)Lf{x(')Ko 

This completes the proof. 

7 Eternal solutions 

Definition 7.1. A weak solution F for min-driven clustering is an eternal 
solution if it is defined on the maximal interval of existence (0, 00). 

Our understanding of eternal solutions is closely connected to the ques- 
tion of how they emerge from clusters of infinitesimal size. Prom the solution 
formula (2.30) we see that as t approaches zero, 



hence Ft{q) — > 1, and this means the relative size distribution Ft always 
converges to a Dirac delta at zero size. A different scaling is needed to 
distinguish solutions through limits as t | 0. 
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7.1 The class of g-measures 

What we will show is that the class of eternal solutions is in one-to-one cor- 
respondence with a suitable space of measures that can be loosely thought 
of as 'rescaled initial data' at t = 0. This correspondence parallels the 
classical probabilistic characterization of infinitely divisible laws via a Levy- 
Khintchine formula. In probability theory, infinitely divisible distributions 
are parametrized by the Levy-Khintchine representation theorem, which ex- 
presses the log of the characteristic function (Fourier transform) in terms 
of a measure that satisfies certain finiteness conditions. In particular [7, 
XIII. 7], a function uj{q) is the Laplace transform L e~'^^F{dx) of an in- 
finitely divisible probability measure F supported on [0, cxd) if and only if 
uj{q) = exp{—rj{q)) where the Laplace exponent rj admits the representation 

V{q) = I ^—^^Gidx) (7.1) 

J[0,oo) X 

for some measure G on [0, oo) that satisfies 

/ (1 A y-^)G(dy) < oo. (7.2) 

J[0,oo) 

(Here a Ab = min(a, b).) As in [11] we call such measures g-measures (short 
for "generating measures," a term motivated by their connection with gen- 
erators of convolution semigroups in probability [7, XIII.9(a)]). Some basic 
analytic facts about Laplace exponents and g-measures are collected in [11, 
Sec. 3]. 

Definition 7.2. A measure G on [0,oo) is a g-measure if (7.2) holds. In 
addition, we say that a g-measure G is divergent if 

G(0) > or / y'^Gidy) = oo. (7.3) 

J[0,oo) 

Here recall we use the notation G{x) = J\Qx]G{dy). The space of g- 
measures has a natural weak topology which is fundamental in our study of 
scaling dynamics. 

Definition 7.3. A sequence of g-measures G^"'' converges to a g-measure G 
as n ^ CO, if at every point x € (0, oo) of continuity of G we have 

G^"'\x)^G{x) and / y'^G^'^Xdy) ^ f y-^G{dy). (7.4) 

J[x,od) J[x,oo) 
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7.2 A Levy-Khintchine formula 

Our analysis of eternal solutions is motivated by a classification theorem of 
Bertoin for Smoluchowski's coagulation equation with additive kernel [1]. 
Here eternal solutions were shown to be in correspondence with divergent 
g-measures. This theorem was generalized to other solvable kernels in [11], 
based on the observation that there is a natural Laplace exponent T]t asso- 
ciated to every solution. For the model now under study, it is determined 
by the g-measure 

xFtjdx) 

and the associated Laplace exponent is 

Vtil) = ^;^^^ = I -^G,{dx). (7.6) 

(Recall that k* from (5.6) is the de Bruijn conjugate of k from (3.13), and 
that K*(0+) = oo if "^pkk log k = oo, k*(0^) < oo if "^p^k log k < oo.) 

Theorem 7.1. (a) Let F be an eternal solution of (4-5). Then there is a 
divergent g-measure H such that Gt converges to H as t I 0. 

(h) Conversely, for every divergent g-measure H , there is a unique eternal 
solution F of (4-5) such that Gt converges to H as t 10. 

(c) The Laplace exponent of H is related to the min history l{t) by 

log^.(g) = f (l - ^-''^'^) V - f^~''^'^^- (7-7) 

To fix ideas, it may help to note that the self-similar solutions are gen- 
erated by the power-law Laplace exponents 

V*{q) = e'\', ^€(0,1], (7.8) 

where 7 is the Euler-Mascheroni constant, which satisfies [14] 

7= / — ds- / ^ds. (7.9) 

Jo -5 Ji s 

This normalization is chosen so the min 

le{t) = t'/\ ^G(0,1]. (7.10) 

The corresponding divergent g-measures are given by 

^^(^^ = f^rrf''' ' ^e(o,i]. (7.11) 
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Proof of Theorem 7.1. In all that follows, q > is fixed, and we use the 
equivalence between convergence of g-measures and pointwise convergence 
of Laplace exponents, as established in [11, Sec. 3], for example. 

1. First, assume that F is an eternal solution. We claim that for all 
q > 0, rit{q) — > ?7^,(g) as t | 0, where tj^, is given by (7.7) and satisfies 
?7^,(oo) = oo. By [11, Sec. 3] it follows that rj^ is the Laplace exponent of a 
divergent g-measure H and Gt — > i/ as t J, 0. 

2. We first recall from Theorem 4.3 that the min history satisfies (4.4) 
whenever < to < *• Regarding t as fixed and to variable, we conclude that 

< l{s) < Cs(i°g2)/Qi ^ < s < 1. (7.12) 

It follows from this, the estimate 1 — e^''-^^'^' < ql{s), and (4.4) that the 
integrals in (7.7) converge, so that ??*(g) is finite for < q < oo. Moreover, 
rj>t:{0) = and ?7*(oo) = Jq ds/s = oo. 

3. We let w = Wt{q) = 1 — Ft{q) = tK'^{t)rjt{q) and use the solution 
formula (2.30) together with (3.13) to write 

log{wK{w)) = -^{Ft{q)) = - / e-«'(^)-. (7.13) 

Jt s 

Adding — log t = J^ ds/s to both sides we find 

(l - e-«'Wj — - / e-«'W^ = logr/,(g) + o(l) 

(7.14) 
as t J, 0. Hence wk{w) ^ tri^:{q), and asymptotic inversion yields 

w ~ tr]4q)K*{tr]^{q)) ~ tr],{q)K*{t) (7.15) 

since k"^ is slowly varying. But immediately this yields r]t{q) -^ 'n*{l) ^-s 
t J, 0, and this finishes the proof of (a). 

4-. We now establish the converse. Let H he a. divergent g-measure 
with Laplace exponent 77* (not known at first to satisfy (7.7)). We will first 
establish that the Levy-Khintchine formula (7.7) defines an appropriate min 
history /, then verify that / defines an eternal solution. 

First, we remark that the definition of the trace admits a natural modi- 
fication for eternal solutions. Given the min history of an eternal solution, 
we define the trace through (3.3). Conversely, we say ^ is a maximal trace 
if A is a distribution function on (0, 00) such that (i) lirriT-^o ^(''") = —00, 
(ii) limT-^oo ^ij) = 00, and (iii) Jq^° TA{dT) < cx) for some tq > 0. In this 
case, the min-history is given by (3.4) with to = 0. 
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Since ??* is the Laplace exponent of a g-measure, 7]^ and 1/ry* are com- 
pletely monotone functions. Thus, there is a positive measure A such that 

(log r?,)' ='k= r e-«V^(dr). (7.16) 

V* Jo 

For every tq > 0, the measure A satisfies the finiteness conditions 

ao^ny>l «"'"C^^(^-)' ^7 17) 

(logr/,) >| ^^^j-g-gr^(^^). (7.17) 

Therefore, we may integrate (7.16) between q and qi £ (0, oo) and rearrange 
terms to obtain 

/•TO fCO 

^ogriM- {l-e~^^)A{dT)+ e-«"^(dT) (7.18) 

Jo Jro 

/•TO /•OO 

= log r/,(gi) - / (1 - e-'?^^)A(cir) + / e-'?iM(dT) := C(gi, tq). 

JO Jto 

The right hand side is independent of q. We let g ^ oo on the left hand 
side, and use 77(00) = 00 to see that /J^" A^dr) = 00. Similarly, we let g ^ 
and use 77(0) = to see that J A{dT) = 00. Thus, A defines a maximal 
trace. Let / denote the associated min history given by / = exp(A^). 

5. Since tq > is arbitrary, we may suppose tq is in the range of / and 
/(to) = Tq. We change variables in (7.18) to obtain 

\ogr,M- f (l-e-^'W) ^+ re-«'W^ = C + logto. (7.19) 

In order to obtain the Levy-Khintchine formula in the form (7.7), we just 
replace / by the rescaling l{s) = l{as) where logo = C + log to- 

6. It remains to check that the solution formula (2.30) defines a solution 
for every t > 0. This is proven by an approximation argument. We consider 
a sequence of finite g-measures G„ that converge to the divergent g-measure 
H. We may suppose that Cn ■= J^ x~^Gn{dx) > n. Let t„ = 1/cn so that 
< t„ < n~^. We will show that the sequence of solutions F^"-* defined for 
t > tn, with initial data given by the probability measures 

Fn{dx) := tnK*{tn)x^^Gn{dx), 

converges to a solution F satisfying (2.30). 

7. Let rjn be the Laplace exponent of G„; then r]n{q) -^ i]*iQ) ^or q > 0. 
Define 

Wniq) ■■= 1 - Fn{q) = tnK*{tn)Vn{q)- (7.20) 
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This is related to the min history l^ determined from F^ by the solution 
formula as in (7.13), namely, 

^ ds / _„j t^\ as 



log{WnK{Wn)/tn)= j (l - e'^'-W) ^ - y g-^'-W^. (7.21) 

Since Wn{q) ^ tnr]n{(l)i<^{tnf]n{(l)) as n — > oo, asymptotic inversion yields 
Wnn{wn) ~ tnTinil), and then (7.21) yields 

logr?.(^)= hm /Vl-e-«'^"'(^))^- re-'^'*"'W-. (7.22) 

Convergence of completely monotone functions also implies convergence of 
all derivatives. Thus, 

It follows that the inverse functions /("' —^ I at all points of continuity. 
Therefore, we may let n ^ oo in the solution formula 



i^)(„\ _,,-! / /°°, -<?/(") (s)^ 



to see that F defines an eternal solution. D 
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